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Abstract 

The equilibrium density profiles in a classical multicomponent plas- 
ma near a hard wall made with a dielectric material characterized by 
a relative dielectric constant e w are studied from the first Born-Green- 
Yvon (bgy ) equation combined with Poisson equation in a regime 
where Coulomb coupling is weak inside the fluid. In order to pre- 
vent the collapse between charges with opposite signs or between each 
charge and its dielectric image inside the wall when e w > 1, hard-core 
repulsions are added to the Coulomb pair interaction. The charge- 
image interaction cannot be treated perturbatively and the density 
profiles vary very fast in the vicinity of the wall when e w ^ 1. The for- 
mal solution of the associated inhomogeneous Debye-Hiickel equations 
will be given in Paper II, together with a systematic fugacity expan- 
sion which allows to retrieve the results obtained from the truncated 
BGY hierarchy. In the present paper the exact density profiles are cal- 
culated analytically up to first order in the coupling parameter. The 
expressions show the interplay between three effects : the geometric 
repulsion from the impenetrable wall; the electrostatic effective attrac- 
tion (e w > 1) or repulsion (e w < 1) due to its dielectric response; and 
the Coulomb interaction between each charge and the potential drop 
created by the electric layer which appears as soon as the system is 
not symmetric. We exhibit how the charge density profile evolves be- 
tween a structure with two oppositely-charged layers and a three-layer 
organization when e w varies. (The case of two ideally conducting walls 
will be displayed elsewhere). 

KEYWORDS : Coulomb interactions, dielectric wall, bgy equa- 
tion, inhomogeneous Debye-Hiickel equation, electric layer. 
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1 Introduction 



The present paper provides new exact analytical perturbative results for the 
density profiles of a classical Coulomb plasma in the vicinity of a polarizable 
boundary. We consider a multicomponent plasma, namely a system made 
of at least two species of moving charges with opposite signs. The linear 
electrostatic response of the wall is described at a macroscopic level by a 
relative dielectric constant e w . (e w is the ratio of the dielectric constants in 
the wall and in the half-space occupied by the Coulomb fluid). The density 
profiles are obtained in a high-temperature (or low-density) limit which is 
realized for instance in an electrolyte solution. 

As shown in Paper II - published just after the present paper - this limit 
is the first-order result in a systematic expansion in powers of the Coulomb 
coupling parameter. This limit can be retrieved from a mean-field approx- 
imation for the first Born-Green- Yvon (bgy ) equation which leads to the 
resolution of inhomogeneous Debye-Huckel equations. For the sake of ped- 
agogy, the present paper is devoted to the mean-field interpretation, which 
should be more familiar to readers interested in chemical physics, and to 
the discussion of the properties of the electric layer. The exact derivation is 
postponed to Paper II where we present two points: first, systematic resum- 
mations of Coulomb divergencies in the framework of the grand-canonical 
ensemble; second, the resolution of the inhomogeneous Debye-Hiickel equa- 
tions obeyed by the auxiliary effective potentials which arise from the latter 
resummations. 

The Coulomb pair interaction u(r;r') between two unit charges located 
respectively at r and r' near the dielectric wall, namely the solution of 
Poisson equation 



where A e i = (e w — l)/(e w + 1) and r'* is the image of r' with respect to the 
plane interface Q. When e w varies from to +oo, A c \ ranges from —1 to 
1. In a dielectric material — 1 < A e i < 1. If the Coulomb fluid mimics an 
electrolyte in a solvent described as a rigid continuum medium, e w is the 
relative dielectric constant of the wall with respect to the solvent dielectric 
constant e s and the interaction potential v w in (|l.2| ) is to be multiplied by 
l/e s . The potential ( |1.2|) may be seen as the sum of two contributions. The 
vacuum or "bulk" potential 





with the adequate electrostatic boundary conditions, reads 




(1.2) 




(1.3) 
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is the solution of (1.1) far away from any boundary or in the vicinity of a 
wall with no electrostatic response (e w = 1). The second term in ( |1.2| ) is 
the interaction with an "image" charge; the latter describes the interaction 
with the polarization charge generated in the material by plasma charges. 
The corresponding self-energy of a charge e a at point r - namely the work 
needed to take one charge e a from the bulk to point r in the vicinity of the 
wall - is equal to e„V^ elf (r) with 

VJ f (r) = -[v w -v B ](T,r) (1.4) 

In ( |1.4| ) the factor 1/2 in the interaction between a charge and its image 
comes from the proportionality between the two charges. In the follow- 
ing, the interface is perpendicular to the x-axis and located at x = 0, and, 
according to ( |1.2j ) , 

V BOli (x) = -A e i^ (1.5) 

When A e i > a hard-core repulsion from the wall must be introduced in 
order to prevent the collapse of each charge with its image. For the sake 
of simplicity, the range b of the repulsion from the wall is chosen to be the 
same for all species in the present paper. Even in the bulk, a short-distance 
cut-off must be introduced in order to prevent the collapse of the system due 
to the attraction between charges with opposite signs. However this second 
cut-off proves not to arise in the densities at leading order in the Coulomb 
coupling parameter inside the fluid. (Indeed, in the first BGY hierarchy 
the variation of the density of every species depends on correlations only 
through an integral and the value of the latter integral at leading order 
in the Coulomb coupling parameter is determined only by the behavior of 
correlations at distances far larger than the short-distance cut-off.) 

For a long time, the short-distance singularity of the charge-image in- 
teraction has prevented one from getting exact results in the case e w ^ 1 
at any distance from the wall for either a generic multicomponent or a 
One-Component Plasma (OCP), namely a system made of only one mov- 
ing charge species in a rigid neutralizing background. The self-consistent 
method introduced by Guernsey || for a plain wall (e w = 1) was general- 
ized for the first time to a case where e w 1 by Alastuey [fflj. This au- 
thor dealt with the OCP near a wall with a repulsive electrostatic response 
(e w < 1) in the weak-coupling limit. In this case the density vanishes on 
the wall and drastically varies over the closest approach distance (3e 2 . The 
mean-field electrostatic potential $(x) created by the charge density profile 
is solution of an inhomogeneous Debye-Hiickel equation where the inverse 
Debye length depends on the distance from the wall and rapidly varies in 
its vicinity. Alastuey solved the equation for the mean- field value of <&(x) 
and produced the corresponding profile density but only for distances larger 
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than the closest approach distance. (For these distances a linearization may 
be performed and the equation for $(x) is a second-order linear differential 
equation with constant parameters). The case e w > 1, where the attractive 
response of the wall makes the density diverge exponentially fast on the wall 
in the absence of any hard-core repulsion, was left unsolved at any distance. 

In Section 2 we introduce a self-consistent scheme for the determination 
of density profiles in a multicomponent plasma from the first BGY equation 
combined with Poisson equation. By using the results of Paper II about the 
solutions of the corresponding inhomogeneous Debye-Hiickel equations, we 
give their formal expressions at any distance from the wall at first order in 
the coupling parameter e D inside the fluid, 

£ d = 7}Pe 2 K D (1.6) 

In (|l.6| ) j3 is the inverse temperature, [5 = l/k B T where k B is Boltzmann 
constant, e is the typical charge in the plasma and n D is the inverse Debye 
length 



k d = knP^e&g (1.7) 

y a 

where is the bulk value of the density for species with index a (and e a 
is the charge of species a). The sum over a runs from 1 to the number of 
species n s . Every density profile takes the form 

p a {x) = pi 9{x - b) e -Wi f (-) [i _ fiec&ix)} (1.8) 

where V^ f (x) is a screened self-energy and <3?(x) is the electrostatic potential 
created by the charge density profile ^ 7 e 7 p 7 (x). (<£(x) is set to in the 
bulk.) The Heaviside function 9{u) - with 6{u) = 1 if u > and 6{u) = 
if u < - describes the geometric constraint enforced by the impenetrable 
wall. The analytic expressions are obtained in Section 3. $(x), given in 
( 3.1?1 ), decays as exp(— k d x) at large distances. V s 3 ^ t (x) may be written as 



the sum 



( x ) = -y l ( k dx; K D b, A e i) - exp(-2/t D x) (1.9) 



V*£f{x) falls off as exp(— 2k, d x)/Ax when x goes to infinity for all values 
of A e i. {k d /2)L given in ( p.9[ ) and ( p. 11 ) arises mainly from the geometric 



deformation of the screening cloud around a charge in the vicinity of the wall 
and remains finite at any distance. On the contrary, — e^(A e i/4x) exp(— 2k d x) 
is the part of the screened self-energy originating from the bare self-energy 
e^Vseii ( |1-5| ) due to the dielectric response of the wall. The second term in 
V s s e ° { (x) was derived for the first time in the case e w < 1 from a phenomeno- 
logical mean-field argument by Onsager and Samaras in 1934 fl"lf . Its con- 
tribution to the density profile is crucial at short distances. When e w < 1 
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(A c i < 0) all charges are electrostatically repelled by the wall, the short- 
distance repulsion range b can be set to zero and the profile density vanishes 
exponentially fast at the contact (x = 0) with the wall. On the contrary, 
when e w > 1 (A e i > 0), all charges are attracted by the wall, b must be kept 
finite and the contact value p a (b) increases as exp [A e i/3e^ exp(— 2k d 6)/(46)] 
when b becomes small. In Section 3 we also derive the profile density in a 
OCP and we compare our result with that of Ref. |fl]. 

Section 4 is devoted to generic global properties of the plasma at the 
interface. In Section 5 we study the case of a plain hard wall (e w = 1). 
The analytic expressions are rather simple and we can investigate the only 
two effects which interplay : the geometric repulsion from the wall and 
the interaction with the electrostatic potential drop <&{x) created by the 
electric layer itself. In the case of a symmetric two-component plasma, we 
retrieve the results of || . In Section 6 the generic properties of the density 
profiles when e w ^ 1 are interpreted in terms of the competition between 
three effects : the two ones already at stake in the vicinity of a plain hard 
wall plus the electrostatic (repulsive or attractive) interaction due to the 
dielectric response of the wall. In particular, we exhibit how the structure 
of the charge density profile evolves from a double layer into a threefold 
layer and then into an inversed double layer when e w increases from the 
value e w = 1. 



2 Self-consistent scheme in the weak-coupling re- 
gime 

2.1 Exact first BGY equation 

The exact density profile p a (x) is related to the Ursell function h a ^ between 
species a and 7 through the first equation of the BGY hierarchy equation, 

d d 

— {\n Pa (x)) = -p— (e a $(x) + e 2 a V scll (x)) 

- f3e a J dr' ( ^ e 7 p 7 (x')h ay (r; r) j ^-(r';r) (2.1) 



In ( |2,1| ) V^ c if(x) is the self-energy (1.4) due to the dielectric response of the 
wall, while <f>(x) is the electrostatic potential created by the charge density 
profile Yl-y 

(x). 3>(x) obeys Poisson equation 

A$(x) = -4vr e 7 p 7 (x) (2.2) 



5 



is uniform in the bulk, since a fundamental property of Coulomb systems 
is the local neutrality relation obeyed by the bulk densities 



(2.3) 



for any value of the Coulomb coupling parameter. Thus, if we redefine 3>(x) 
as the difference between the electrostatic potential created by the charge 
density and its bulk value, tends to zero when x goes to +00. 

We recall that, in the case e w = 1, where there is no image forces, the 
density is merely uniform in the zero-coupling limit. In a plasma with no 
charge symmetry the potential drop <3?(x), which is determined from the 



charge density profile Y^-y e iPi{ x ) through Poisson equation (2.2), does not 
vanish. Moreover, in the weak coupling regime, it is of the same order as 
the pair-correlation contribution to dp a /dx in the BGY equation, as shown 
a posteriori by our explicit calculations displayed in Section 3. By a mean- 
field scheme Guernsey || closed the second BGY equation in a weak-coupling 
limit and found that the zeroth-order pair correlation is calculated with uni- 
form densities in a semi-infinite space. Thus he obtained coupled equations 
for dpa/dx and <I>(x) and calculated J2 a e a p a (x) as a double integral. (The 
resolution was not performed for every p a {x) in the case <&(cc) 7^ though 
it might have been done.) The density profiles p a (x) were studied only in 
the case of a symmetric two-component plasma. The first-order correction 
to their bulk value in the weak-coupling regime was calculated by Jancovici 
H as follows. Because of the charge symmetry specific to this system, the 
charge density profile and subsequently the electrostatic potential difference 
with the bulk $(x) vanish at any distance from the wall. Then the gradi- 
ent of the density p a {x) of species with charge e a given by the first BGY 



equation (2.1) is determined at leading order only by the zeroth-order pair 
correlation ; the latter is calculated in a mean-field approximation for the 
direct correlation function with the same result as that found by Guernsey 
1- 

When e w 7^ 1 the methods introduced in the case e w = 1 cannot be 
generalized straighforwardly, because the fast variation of the density in the 
vicinity of the wall prevents one from using mere linearizations. In order to 
circumvert this difficulty we introduce the following scheme. 



2.2 Mean-field Ursell function 

In the first BGY equation ( ^l| ) ]P ^jPj(x')h a -y(r; r') is the excess charge den- 
sity of the screening cloud around a charge e a located at r, the excess charge 
being calculated with respect to the charge density profile £V e 7 p 7 (x). The 
electrostatic potential created at r" by the charge e a and its screening cloud 
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IS 



^ > cxc,a(r;r // ) = e a v w (r;r") + J dr' e 1 p 1 (x')h ai (r; r')j v w (r';r") 

(2.4) 

A mean-field approximation amounts to assuming that 

«, Q (r;r")=e^ MF (r;r") (2.5a) 
h™ F (r; r') = -(3e a e 7 <f> MF '(r; r') (2.5b) 

(At leading order in the parameter e D , the long-range Coulomb interaction 
prevails and the short-distance repulsion between particles is not involved 



in (2.5)). By inserting these approximations into the definition ( [2.4D , we 



obtain the well-known mean-field equation 



(r;r") = « w (r;r")-/3 J dr' e?p 7 (x')^w(r; r')cf> MF (r'; r") (2.6) 



Then the mean-field approximation of the integral in (|2.1| ) proves to be 
equal to 

-^[</> MF --w](r;r) (2.7) 

Indeed, the integral in (^]l]) can be rewritten by means of the trick involving 
the Dirac distribution <5(r — r"), 

J dr'g(r;r')^(r';r) = J dr" 5(r - r")-^ (J dr' g(v"- r')/(r'; r)) (2.8) 

Moreover, according to fll.ip and ( |2.6| ), MF (r;r') is the Green function 
of the operator A r / — 4-7T/5 ^ 7 e^p 7 (x'). Since the latter operator is self- 
adjoint, the real function (j) MF (r; r'), as well as v w (r; r'), is symmetric under 
exchange of r and r' when r and r' are in the same region. Then ( |2.6| ) can 
be used. Since a symmetric function h(r;r') = h(r';r) obeys the identity 
d[h(r;r')]/dx\ I=T > = (l/2)d[h(r;r)]/dx, we get the result Q. 



Finally, according to the definition ([LJ) of V seU (x) and since p a (x) tends 
to p„ when x goes to +oo, the mean-field density profile p^ F (x) proves to 
read 

p™ F (x) = 9{x - b)pl exp [-(3elV;: it {x) - 0e a <S> MF (x)] (2.9) 

with 

KS(x) = J(^"-«B)(r,r)- lim I(^'-^)( r ,r) (2.10) 
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Meanwhile the coupled equation for the electrostatic potential <fr MF (x) is 
(2.2) where p 1 {x) is replaced by p* IF (x). In ( |2,9[) the argument in the 



exponential may be interpreted as (3 times the work given by an operator to 
the system in order to put a charge e a into the Coulomb fluid at r, make it 
cross the potential drop $ MF (x) from x to +00 and then get it back from 
the bulk. In the following, we set 

z™ F {x) = 9(x - b)pl exp [-(3e 2 a V;: it (x)\ (2.11) 
2.3 Linearization of the $(x) contribution 

Explicit calculations can be performed if the contribution from <£ MF (x) to 

p™ F ( x ) = z ^ F {x)e-^ MF ^ (2.12) 

is linearized, 



f%*<°*{x) = z A a IF (x) [1 - 0e a $ MF ^(x)\ (2.13) 

Such a linearization is allowed only if $ MF (x) does not become infinite in the 
vicinity of the wall. This is indeed the case, because $ MF (x) is created by 
the charge density and the latter has no singularity thanks to the hard-core 
repulsion from the wall. The absence of divergency in <3?(x) near the wall is 
also checked in the systematic approach of Paper II. In the following we will 
show that the screened self-energy V^ f (x) diverges when x goes to zero and 
no linearization can be performed for it. 

By inserting the linearized mean-field expression for the densities into 



Poisson equation (2.2) we find that 

A r - AttPY,^ z j F ( x ) $ MFM (x) = -4vr^e 7 4 /F (x) (2.14) 



( 2.14| ) can be viewed as some kind of partially linearized Poisson-Boltzmann 



equation in an inhomogeneous case. As a consequence of ( |2.14 ), 



<S,m F m (x) = f dr > «^V, r) £ e 7 < F (*') (2.15) 
J 7 

where <j>^ F ' (r, r') is a Green function solution of 

[A r - 4 (1 + U(r))] r' Un (r,r') = -4n5 (r - r') (2.16) 

with 

[/(r)4E e ?[< F W-^] ( 2 - 17 ) 

K D 7 

0AfF,hn j s ^ e solution of ( p,16|) which satisfies the same boundary conditions 
as v w (r, r'). 
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2.4 Solution of the inhomogeneous Debye-Hiickel equation 
at leading order 

A formal solution of the inhomogeneous Debye-Hiickel equation ( p. 16 ) is 
given in Paper II. An e^-expansion is devised and we show that 

^ 1FAia (r,r>) = K JW( KD r, KD r') [1 + O cxp (e D )] (2.18) 



In ( 2.18] ) O cxp (e D ) denotes a function of order e D - possibly multiplied by 



some power of lne D - which decays exponentially fast at large distances 
over a scale k^ 1 and which remains bounded by a function of fie 2 /b for all x 
larger than the closest approach distance b to the wall. k d 4>^ is the solution 
of the homogeneous Debye-Hiickel equation 

[A r - 4] k d ^ (k d t, k d t') = -A<k5(v - r') (2.19) 

with the same electrostatic boundary conditions as v w (r, r') : 4>(°> is contin- 
uous in all space and 

50(o) q7{o) 

lim e w — — (k d t, k d t') = lim — — (k d t, K D r) (2.20) 
cc->o- ox x^0+ ox 



Equation ( 2.18; ) also holds for <f> MF defined in ( |2.6| ) and which obeys equation 
( 2.16|) where z^ IF {x) in U(r) is replaced by p-y(x). Thus the screened self- 



energy defined in ( |2.10| ) is also determined at leading order in e D 

Eventually, the density profile at leading order in e D is given by ( |2.11 ) 
and ( 2.13; ) with 



<t> {0) - 0b| (« D r, k d t) (2.21) 



and 



^,iin (x) = J dr ' KD ^)^ DriKDr ')Y^e^ F {x') (2.22) 



In (]2.21 ) cp B denotes the solution of the homogeneous Debye-Hiickel equation 
in the bulk, namely the solution which vanishes when |r— r'| becomes infinite, 
as well as v B (r, r') defined in (|1.3|). 



3 Explicit density profiles 

3.1 Solution of the homogeneous Debye equation 



Equation ( pltjp can be solved because it is changed into a second-order 



differential equation by taking the Fourier transform in the direction parallel 
to the wall, 

^ (0) (x,x',q)= / dye i(i -yft 0) (x,x',y) (3.1) 
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In (3.1) we have used the dimensionless variables and y — K D y, 

where y is the projection of r onto a plane parallel to the wall. As recalled in 
Paper II, since K D (p^ obeys the same boundary conditions as v w , for x > b 
and x' > b 



(/> (0) (x,x',q;6, A c i) = <f> B (\x - x'|,q) 

+ Z(q; b, A cl )h+ w {x + x' - 2b; q) (3.2) 



with 



<f) B (\x - x'\,q) 



2tt 



VT+q 2 



-\x-x'\yjl+q 2 



(3.3) 



ht w (x + ~x> - 26; q) = V^-H 

VI + T V 1 + Q + nl 



and 



£(q; b, A e i) 



1 - A el e- 2 ^ 



Vi + g 2 + |q| 



1 - A c ie- 2 i b 



(3.4) 



(3.5) 



3.2 Screened self-energy 



The screened self-energy ( 2.21 ) can be written as 



0e 2 a V^ ( (x)=e a L(x-b;b,A e 



with e a = (l/2)/?e 2 k d and 



L(x - b;b,A c \) = 



d 2 q r 



(2^ 



According to (|3^ 

L(u;b, A el ) 



(i 2 q 

(2vr) 2 



Z(q;6, A e i) x h+ w (2it;q) 



(3.6) 



(o) (x, 5, q; 6, A e i) - <£b(£, 5, q) (3.7) 



(3.8) 



where Ahw is given in ([O]) . By using the change of variable t = \J\ + q 2 
we get 



L(x — b; b, A e 



+oo 1 — A e i 
dt 



-2{x-b)t 



t + - Aeie-^v^ 



(3.9) 
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The successive changes of variables t = t' + 1 then t' = vj 2(x — b) 
allow to show that 



L(x-b:b, Ani) ~ 



-2(x-b) 



x^+oo 2(x — b) 



(3.10) 



when 



If b 7^ the integrand in ( |3,9| ) behaves as 1/i 2 times exp —2{x — b)t 

t goes to +oo and L(x — b, b, A e i) is finite for all values of x even when x 
approaches b. If b = the integrand vanishes as exp(— 2xt)/t 2 for large i 
when A e i = but it behaves as — A e i exp(— 2xt) if A c i ^ 0. Subsequently 
for 6 = 0, the integral diverges at x, = b = when A e i ^ 0. By subtract- 
ing the dangerous asymptotic behaviour I^t) = — A c iexp [— 2xt] from the 
integrand of L and by performing f^~°° (it I as (i) we get 



-21- 



L{x - b; b, A e i) = -Agi-— + L{x; b, A ei ) 



2x 



where L(x; b, A e j) remains finite even when x = b = 0. 



(3.11) 



3.3 Electrostatic potential drop 

In order to calculate the electrostatic potential drop ( |2.22| ), we notice that 
(x, x', q = 0; b) = 2vr [e - '* - *'! + e ^ + ' x '- 2 ^} (3.12) 
and we rewrite z^ F (x) as 

Za F {x) = Pi [1 + ^o(x; e a , A c i)] [l - £ a £(x; &, Aei)] (3.13) 

with 



w (x;e a , Ad) = exp 



4x 



1 



(3.14) 



According to the bulk neutrality relation ( |2.3| ) , the constant term in z^ F {x) 
gives a vanishing contribution to fl2.22|) . wq(x; e a , A e i) proves to contribute 
from order e D lne D to the integral in ( |2.22j ). Indeed, let us consider a func- 
tion f(x';x,b) which is bounded for all x! > and b > 0. If b > k^ 1 , then 
for all x > b Pe^/x < K D f3e 2 Y and 



+oo 



dx' w (x';e a , A e i)/(x'; x, b) 
dx' 



+°° w A e i/3e 2 



4s' 



°e- 28 '/(5 , ;x,6) 



(3.15) 
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In the case b <C k^ 1 , let us introduce the length I such that /?e^ <C I <C k^ 1 . 
For all x in the range b < x < I, n D x <C 1, while, for all x > I, /3e 2 /x <C 1. 
Then at leading order in e D 

k d / dx' w (x;£ a ,A ei )f(x;x,b) 
Jb 



~ lim lim 



dx' 



cxp 



4x' 



- 1 



/(x' = 0;x,6 = 0) 



+k d 



45' 



e" 2 * f(S>;x,b = 0) 



1 + 0(6) 



(3.16) 

where C(6) denotes a term of order 6. After the change of variable x 1 = 
x^/Se 2 , the first integral in ( p.!6| ) proves to be of order e a , as well as the 
second integral. Both integrals have a logarithmic dependence upon I and 
the respective ln(Z//3e 2 ) and Iu(k d I) terms combine so that the sum of the 
two integrals starts at order e D In e D . 

Eventually, we get for $ MF,lin (x), denoted by $(x) in the following, 

= P^M^nnx; e 7 , K D b, A el ) x [1 + O cxp {e D )} (3.17) 



where ©exp^o) is defined after (|2.18[) . In ( |3.17D M 7 = M + [M 7 - M] with 

i r+oo 



M(x;b, A e 



2 A 



e -|2-u'[ + e -(S+u'-»)j £( u /.g jA 



(3.18) 



and M 7 — M is the e D -expansion at orders lne 7 and (e 7 )° of the integral 
1 1 f + °° 



2 £ 7 Jl 
3.4 Density profile 



e -\x-u'\ +e -(x+u'-2b) 



w (u , ;£ 7 ,A el/ 



(3.19) 



The previous expressions are inserted in ( |2.13| ) with the result 



p Q (x) =p B J(x - b) exp ( A el ^-e~ 2 ^ x 



X<1 



1 



e 2 a L (k d x; K D b, A el ) 



47T/3 



e a 2~ / PyG.^M^^K D x; £ 7 , A e 



,2 Z_f r 7 7 



(3.20) 
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where e 1 = (l/2)Pe*K D . M 7 is defined in ( gig ) and flSlSD , while L is given 



by ( |3.9| ) and (|3.11| ). L(x; 6, A e i) decreases exponentially fast as exp(— 2x)/2x 
when x goes to oo, while M 7 decays only as exp(— x)/x. 
We give more explicit formulae in the regime 

r] = K D b<l (3.21) 



whatever the value of j3e 2 jb may be. When e D <C 1, according to (|1.6|) 
and ( |1.7|) , the mean interparticle distance a is smaller than k^, 1 , a < k^ 1 
and the condition ( |3.21 ) will be fulfilled if b <C a - for instance if b is of 



the same magnitude as the hard-core diameter of charges which itself is far 
smaller than a. L(x;n D b, A e i) is bounded for every x, even when r/ = K D b 
vanishes, and it can be expanded in powers of r\. According to (|3.9| ) and 
( 3. lip , L(x;r] = 0, A e i) is directly given by 



L(x; ry = 0, A e i) = / dte" 2 ^' clJ (3.22) 



On the other hand 



M 7 (x; e 7 , r; = 0, A e 



M(x;r7 = 0, Ad) + (lne 7 )A el -e- £ - A cl / 7 (x; A cl , /3^/6) (3.23) 



where 
M(x;r/ = 0,A cl ) = 

while, according to (|3,16| ) 
1 



dt 



-2tx 



2te 



1 - (2t) 2 



1 " A c 2 t 



t + Vt 2 - 1 ) -A c i 



(3.24) 



iJx; A e i,/3e 7 /6) 



x <^ 2 



,4 



/3e 2 A cl 



.,2 

46 



+ m|i^L]+2C-] 



+ 4 6 



-2x 



e x Ei(-x) - e 3i Ei(-3x) 



+ ln3 
(3.25) 



In ( p.25| ) C is the Euler constant, Ei(— u) denotes the Exponential-Integral 
function defined for u > as Ei(— u) = — f^°° dt exp(— t)/t and 



A(u) = - [e u - 1] - Ei (u) 



11 



(3.26) 



A(u) arises from the integration of exp(— /?e 2 A e i/4u') — 1 in ( p. 19 ). 
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3.5 Interpretation : competition between three effects 



The profile density is ruled by the competition between three kinds of effec- 
tive interactions, as exhibited by rewriting the density profile by means of 



D, dpD and (|37Tl|) with the result 

x {l " 0e 2 JjL{K D x; K D b, A e i) - (3e a ${x) + O cxp (e 2 D )} (3.27) 



where O exp (e 2 D ) is defined after ( 2.18| ). The interpretation of ( |3 .27 ) is the 
following. 

First, e a &(x) is the interaction between a charge e a and the charge 
profile density in the electric layer. The other two interactions, which are 
proportional to , are the two parts of the screened self-energy. 

Second, exp(A e i/3e„ exp(— 2k d x)/'1x) is the effective Boltzmann factor 
associated with the part of the screened self-energy created by the elec- 
trostatic response of the wall. The effect of the corresponding attractive 
(A c i > 0) or repulsive (A e i < 0) interaction with the wall cannot be lin- 
earized. Indeed, when the dielectric wall is repulsive, the density vanishes 
as exp(— | A c i|/3e^/4&) when b goes to zero. Since the hard-core repulsion is 
spurious when A e i < 0, we can set b = and 

p a (x = 0, Ad < 0) = (3.28) 

On the contrary when A e i > 0, the density blows up as 

p a (x = b, A cl ) ~ o pi exp l Ael(3e }f KDb \ [1 + 0(e D )} (3.29) 

Third, the part of the screened self-energy which exists even in the ab- 
sence of any electrostatic property of the wall (namely even when A e i = 0) 
is (k d /2)L(k d x; K D b, A c \); it arises from the "geometric" repulsion caused 
by the mere presence of the wall. Indeed, by deforming the screening cloud 
(with a net charge — e a ) surrounding any charge e a , the wall hinders the 
stabilizing effect of Coulomb interactions, as exhibited clearly in subsection 
5.3 about the plain hard wall. 



3.6 Limiting case of the OCP 

The density profile for a one-component plasma (OCP) can be derived from 
the expression obtained for a two-component plasma by the following trick 
already tested in Refs. [§, ||]. In order to describe a OCP where moving 
particles carry a positive charge e in a neutralizing uniform background with 
density p, we start from a two-component plasma with e+ = e and p+ = p 
and we take the limit where e_/e+ vanishes while e-p- = — e+p+ is kept 
fixed. 
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The general expression ( [3.20 ) tends to the limit 
~A el /3e 2 e- 2KDX ~ 



Pocp = pO{x-b) exp 



4.r 



jl - ^/3K D e 2 [L(k d x; K D b, A c i) - M ocp (k d x; (3K D e 2 /2, K D b, A e i)] 



(3.30) 



where k d = ^AirPpe 2 , L is defined by (|3T9|) and ( |3.11[) , while M OCP is equal 
to the function M 7 defined in ( |3.18 ) and ( |3.19| ) with e 7 replaced by K D (3e 2 /2. 
According to ( 3.17 ) , the potential difference with the bulk is 

1 



^ocp(x) 



-K D eM OCP (K D x; (3K D e /2, K D b, A e i) 



We recall that ( |3.30| ) can also be rewritten as 

PocP {x) = P 9(x- b) e-^KSM [i _ (3e<S> OCP {x)] 



(3.31) 



(3.32) 



where V s %(x) is given in (|i~9l) . Moreover, according to ( ^22|) may 

be rewritten in the case of the OCP as 



1 

2.1b 



dx 



e -\x-x'\ + e -(x+x'-2b) 



(3e$ OOP (x;b) + O exp (e 2 D ) (3.33) 



3.7 Comparison with previous results 



The expression ( |3,30| ) for the OCP is valid at any distance from the wall 
and for any value of A e i and n D b. Let us see how it can be compared with 
the expression (3.21) in Ref.(l[ obtained in the case A e i < and 6 = 0, 
and for distances x 3> (3e 2 . (The expressions in Ref.[|J will be denoted by a 
superscript *.) 

Alastuey starts from the BGY equation (|2.1| ) and directly replaces the 
Ursell function by its Debye approximation —l3e 2 K D ^ ' (n D r, k d y'; A c i). The 
equation (3.17) in Ref.[jl| - which is analogous to our equation (|2.14| ) - 
involves the screened self-energy calculated at leading order in e D - as in 
the present paper - and given by 



V;° c *(x) = £ D L(K D x;A e 



(3.34) 



However the equation (3.17) in Ref.|IJ is solved only for distances x 3> 
f3e 2 which are large enough to allow one to replace U(r), defined in our 



equation (|2.17[) , by zero. The corresponding approximated &* OCP is solution 
of 



dx 2 



(3e 



x > /3e 2 



(3.35) 
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while, according to the Poisson equation which relates d 2 <&* OCP / dx 2 and the 
charge density e [p(x) — p], 

p* OCP (x) = p [e-P^n*) - (3e% CP (x)] x » (3e 2 (3.36) 

(( |3,36j )is the formula (3.21) given in Ref. |l[].) According to the result ( |2.18| ) 
of our analysis of the non approximated equation ( 2.16j ) (see Paper II), the 



expression of &% CP , solution of ( |3.35 ) and given in Eq.(3.20) of Ref. 



coincides with our formula ( |3.33| ) valid at any distance. On the other hand, 
the result ( |3.36| ) does coincide with the x 3> j3e 2 limit of the expression 
( 3. 32] ) which is valid for any 



x. 



Our ability to handle with all distances relies on two progresses with 
respect to the approach of Ref.[jl]]. First, we extract from the expression of 
e aV s S ci*( x ) the part which diverges as the bare self-energy A e \e 2 I /(4:x) when 
x approaches the wall and which enforces the vanishing of the density at the 
wall when 6 = 0. Second, we are able to disentangle this short-range effect 
from the long-range exponential screening through our systematic method 
of expansion introduced in Paper II for the solutions of the inhomogeneous 
Debye-Hiickel equations at stake. 



4 Generic global properties 
4.1 Potential drop 

First we recall that in the generic case the local neutrality valid in the bulk 
is destroyed near the wall, 



e a Pa(x) / 



(4.1) 



and there appears an electric layer which is responsible for a potential drop 
&(x) between each point and the bulk. The expression of <&(x) is given 
in ( |3.17| ). M in the decomposition of M 7 is a positive function, whereas 
M 7 — M may have any sign. Thus the profile of $ depends on the composition 
{e 7 , /j 7 } 7=1) _ ns of the Coulomb fluid and on the value of A e i. For instance, 
when K D b -C 1, 



= 0) = -t^L £ p B e 3 (1 _ Ae2l)J(Aclj 
Ad " 



+ 



In (3K D pe 2 1 Aei|/4) - 1 + 2C + A 



4b 



+ 0{e D -K D b)\ (4.2) 
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In (O) 0{e d] K D b) denotes a term of order either e D or K D b, and 



J(Ad) - 4(l + A c , + A;,) {-^ + ' 1 + 2A °'>'" 3 




An example for the profile <&(x) is drawn in Fig. |l|. 

However the local neutrality ^2 a e a p a (x) = holds in the specific case 
of a charge-symmetric plasma in a symmetric state, for any strength of the 
Coulomb coupling in the fluid. We use the following definitions. A charge- 
symmetric fluid contains equal numbers of positively and negatively charged 
species and the set of charges is invariant under inversion of charges. In the 
special case of a two-component plasma made of charges +e and — e, the 
latter charge symmetry combined with the neutrality relation ( |2.3|) implies 
that p\ = p*L- On the contrary, for a charge-symmetric plasma with at least 
four species p& ^ p B _ a in the generic case. However, in some situations (for 
instance when two different salts made with monovalent ions are dissolved 
in water) the system is prepared in a symmetric state; the bulk density 
parameters are chosen to satisfy 

PT™ = P-T (4-4) 

In such a symmetric state, the symmetry of the Hamiltonian under inversion 
of charge signs enforces that at any point x 



_ ,(%'°{x)=0 (4.5) 

a 

and, subsequently 

$ s y m ( x ) = (4.6) 

In a symmetric state, a charge-symmetric Coulomb fluid does not build any 
charge density profile or any electrostatic potential difference with the bulk. 

4.2 Global charge 

A dielectric wall remains globally neutral in the presence of a Coulombic fluid 
and may only acquire macroscopic multipoles depending on the geometry of 
the dielectric sample. As a consequence, we expect pl| that, since e w < +oo, 



o 



+oo 

dx ^^e a p a (x) = (4.7) 



We recall that <fr(x = 0) ^ means that the dielectric layer carries a nonva- 
nishing dipole though its net charge is zero. 
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The global charge of the system calculated at leading order e D indeed 
vanishes in agreement with fl4.7|) . The result is readily obtained by writ- 
ing for p Q (x) the structure ( |2,13| ) combined with ( 2.22| ) and by using the 
following property. According to (|3.15| ) and ( 3.16 ), if f(x) is an integrable 
function which is bounded for every x > 0, 



r+oo r+oo 

/ dxz™ F {x)f{x) = pl dxf{x) X [l + 0(e D \ne D )} (4.8) 
J o Jo 

Eqs. (|2~13|) and (gj) lead to 

r+oo r+oo 

/ dx S2e a p a (x) = dxS2e a Za IF (x) 
Jo Jo 



(4.9) 



<3?(x) is given at leading order e D by ( |2.22| ), and the property 



implies that 



P+OO p 

/ dx / dy K D ^°\K D x,K D x',K D y) = — (4.10) 
Jo J K D 

r+oo a r+oo 

/ dx$(x) = ^- / dx'V e 7 < F (x') (4.11) 
Jo 4io „ 



Combination of ( [4.9j ) and ( 4.11| ) implies that (4.7) is indeed satisfied at 
leading order in e D . 

4.3 Contact theorem 

Finally, we turn to the so-called contact theorem which gives the difference 
between the bulk thermodynamical pressure P B and the kinetic pressure on 
the wall k B T^2 a p a (x = b). As shown for instance in Ref. |J, 

2 



PP B =J2pa(x = b)-27rA el f3 



r+oo 



&aPa (•£) 



-(*) 



— 3 ^ dx j dx j dy — - — — (x, x',y) 



(4.12) 



dx 



X 



^ e a e 1 p a (x)p 1 (x')h ai (x, x', y) 



a,7 
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Since the global charge in the vicinity of a dielectric wall vanishes (see ( [4.7| ) ) , 



the second term in the r.h.s. of ( 4,12 ) is equal to zero 



The contact theorem implies that compensations between the various 
terms in the r.h.s. of fl4.12j ) ensure that the bulk pressure is independent 
from b as well as from e w ; in other words, the bulk pressure is independent 
from the specific forms of the interactions between particles and the wall, 
whether the latter interactions are purely geometric repulsions or coulombic 
couplings. 

The above compensations can be checked at first order in the coupling 
parameter e D , as shown in Appendix. On one hand, up to order e D , the 
bulk pressure P B is just the sum of the ideal-gas pressure plus the Debye 
correction || 

pB = - 1| + ° K lne -) ( 4 - 13 ) 

a 

where p is the order of magnitude of the p^'s. On the other hand, we 
calculate the r.h.s. of fl!.12| ) by using the fact that the density profiles at 
first order in e D obey the first BGY equation ( |2.1| ). In order to handle the 
contributions from the screened self-energy in p a {x) properly, the integrals 
are performed by using properties similar to ( [4.8[ ) and derived from ( |3.15| ) 



and ( 3.16 ). Finally, after compensations of terms involving the dielectric 
response of the wall, the ideal-gas pressure arises and the remaining term is 
reduced to the difference between the kinetic pressure at the contact with 
a plain hard wall (e w = 1) and the ideal-gas pressure. The latter difference 
involves only the explicit value of the screened self-energy at the contact 
x = b with a plain hard wall. This value is independent from b and gives 
the term of order e D in the bulk pressure ( 4.1 3| ) . 



5 Case of a plain hard wall (e w = 1) 
5.1 Explicit formulas 

In this case the profile density is ruled only by the competition between 
Coulomb interactions in the fluid and the geometric deformation of screen- 
ing clouds by the impenetrable wall. According to ( 3.17| ) and ( |3.20|) , the 
expression of the density profile is reduced to 

Pt w (*) = Pi {l - PelfL™ {k d {x - b)) - (3e a ^(x)} (5.1) 



In (|5.1| ), L HW can be explicitly calculated from (O) with the result (see 



V 



in) 



-2u 



11 1 

2u u 2 2v? 



-K 2 (2u) 
u 



(5.2) 
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where K2 is a Bessel function. According to ( |3.19| ), M HW is reduced to 
M HW (x). Thus M HW is independent from the species 7 and, according to 



<S> HW {x) 



(« D (x-6)) 



(5.3) 



where M HW depends only on k d {x — b) since it is defined in terms of 
L HW (k d x; K D b) through ( p8j ) and L HW (k d x; K D b) = L hw (k d (x - b)) ac- 
cording to ( |3.11| ) . Its expression at u = k d (x — b) coincides with ( |3.24j ) when 
A e i is set to zero, 



M HW (u) 



dt 



-2tu 



2te~ u 



The large-distance behaviour of 



M 



(u) ~ — 



1 - (2t) 2 

1 reads 
ln3 - 2 + 



t + Vt 2 



(5.4) 



7T 

V3 



(5.5) 



Since M HW is a positive function, § HW {x) has the same sign at all distances 



from the wall. This sign is determined by ^ Py e ^, 



U™(*)<0 



(5.6) 



5.2 Electric layers 

Near a hard wall the charge density profile takes the simple form 

a 

= -0 (X>" e ') T {Kd[x " 6)) " ^ Mx " (5 - 7) 



where —M HW is proportional to according to (5^). We have checked 
that d5^) agrees with the result (32) in Ref. [§. Since L HW (0) = 1/3 and 
M HW (0) = [In 3 + 1 - tt/V3] /8, L mv (0) > M HW (0), while L w and M HW 
are positive functions of x which decay respectively as exp(— 2k d x)/x and 
exp(— k d x) when x goes to infinity. Therefore the expression ( |5.7|) implies 
that the charge density profile is at least a double layer. 

In a charge-symmetric plasma in a symmetric state, the local neutrality 



(4.5) and the vanishing of <J>(x) ( |4.6[ ), which are valid whatever the strength 
of the coupling inside the plasma may be and for any value of the densities, 
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are retrieved at first order in e D from our expressions. Indeed, according to 
(5.7) and (|5.3|) , Yla e aPa( x ) an< ^ & HW {x) are both proportional to X^P^^ 



and this combination vanishes in any charge-symmetric plasma in a sym- 
metric state. In a plasma with an even or odd number of species, for a 
particular set of densities which satisfies the constraint 



E 



fie* = (5.8) 



the properties Q4.5|) and Q4.6| ) happens to be valid at first order in e D . 

In a plasma which is not in a symmetric state, the charge density profile 
does not vanish and the sign of the charge density at the wall is fixed by the 
sign of ^2 p~e* (in the considered weak-coupling regime), 



Ep" e 7 E e ^( x = 6 ) <0 ( 5 - 9 ) 



The inequality (|)Jj) implies that if the magnitudes of positive charges is far 
larger than those of negative charges, then the layer at the contact with the 
wall is negatively charged. Moreover, the combination of ( [5.6D and ( |5.9| ) 
implies that 

^ HW (x = b)J2e a Pa W (x = b) >0 (5.10) 



5.3 Repulsion from the wall for the total particle density 



According to the bulk local neutrality (2.3), 

.3 
8tt 



E p" w {x) = E & - t LHW - 6 )) ( 5 - n ) 



Since L HW is a positive function, the total particle density is lower than its 
bulk value at any point, 



a a 



The total particle density undergoes a repulsion at any distance. 

The repulsion from the wall also operates for every particle density when 
the electrostatic potential &(x) vanishes, namely in the case of a charge- 
symmetric plasma in a symmetric state or in a plasma where the composition 
happens to satisfy ( |5.8j ). Indeed, according to ([0]), at any point 

P" W (*)k , =0 <pg (5.13) 
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The repulsion from the wall for every species when the potential drop with 
the bulk vanishes is interpreted as follows. According to (|2.9| ), when $(x = 
b) = 0, the ratio p" w { x = &)/Pq is onr y determined by the screened self- 
energy which is the difference between the values of (e^,/2)[0 — v B ](r, r) at 



the wall and in the bulk. According to (5.1), this difference is equal to 
{e 2 a /2)n D L HW (0), which is positive. Thus, when there is no potential drop 
between the wall and the bulk, the immersion free-energy in the bulk is lower 
than its value at the hard wall : the charge surrounded by its screening cloud 
with global charge of opposite sign is more stable in the bulk than at the 
plain wall. In other words, for all species Coulomb screening is less efficient 
when polarization clouds are deformed by the presence of the hard wall. 
As an illustration, we consider a symmetric two-component plasma made of 
charges e and — e. According to charge-symmetry p+(x) = P-{x) and the 
profile density is drawn in Fig. ||:a. 

5.4 Particle density profiles 

The case of plasmas in a charge symmetric state has been discussed in the 
previous subsection. Here we consider the generic case where ^ 7 P^e^ / 0. 
Then, already at leading order in e D , & HW (x) ^ and ^2 a e a p a (x) ^ so 
that the electrostatic potential created by the electric layer interplays with 
the geometric repulsion from the wall. 

If J2 7 Py e ^ > then ® HW (x) < for all x according to (|5^ ). Thus for 
all negatively charged species a~ at any point, according to ( |5.1| ), 

p™(x) < p B a _ (5.14) 

because the geometric and electrostatic effects are both repulsive for them. 
For positively charged particles, since L HW (x) decays faster than & HW (x), 
the attractive effect of the potential drop overcomes the wall repulsion at 
sufficiently large distances, and p"Y{x) — p® becomes positive a priori at 
least at some distance xq before decaying to zero when x goes to oo, 

p H J{x)- p B a+ >U for x>x (5.15) 

The result of the competition between the two effects is given by the density 
on the wall, which reads 



p™(x = b)=plll-±f3K D 



7T 



(5.16) 



The sign of p^ w (x = b) — p^ depends on the particular composition of the 
plasma. 
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More precise results about the layer structure of p" w (x) — p^ 
obtained in two special cases. First, in a two-component plasma made of 
charges e + and e_, the local neutrality in the bulk Q2.3]) enforces that 



Pl w {x)=pl 



L HW (k d (x - b)) 



1 + 



M HW {K D {x-b)) 



(5.17) 



The argument displayed after ( |5.7| ) shows that if e + > |e_|, p+ w (x) < p\ 
in a strip b < x < xq whereas p I f w (x) > p+ for all x > xq. As a conclusion, 
at leading order in e D , the wall repulsion still overcomes the electrostatic 
attraction arising from & HW for the positive charges near the wall and the 
profile density p" w (x) — p\ has the structure of a double layer. This can be 
seen in Fig. §:b. 

At last, we briefly discuss the case of a three-component plasma with 
> and which is made for instance of species e\ > 0, ei > and 



e% < 0. Then p%(x) < pf according to ( 5.14| ). However, according to ( |5.16| ) 



the composition {p®, e 7 } 7= i j 2,3 of the fluid may happen to be such that the 
electrostatic attraction of positive charges (proportional to e a ) overcomes 
the geometric repulsion from the wall at x = b (proportional to e^) for the 
species which carries the positive charge ej with the lowest magnitude, for 
instance &i. Then piix = b) may happen to be larger than pf in spite of 
the wall geometric repulsion. If p2{x = b) > pf > according to (5.9) and the 
neutrality condition (|2.3| ), pi(x = b) < pf , and pi(x) — pf contains at least 
a double layer. 



6 Generic local properties 

6.1 Large-distance behaviours 

As shown in Section 2 the density profile takes the form The total 

screened self-energy e^V s s ^ { (x) is written in (|3.6|). Its decay at large distances 
e 2 a exp(— 2k d x)/4x (given by ( 3.10| )) is independent from A e i and is positive : 



far away from the wall, the screened self-energy is a repulsive effect, even if 
the electrostatic response of the wall upon one charge is attractive (A e j > 0). 
The contribution from the complete screened self-energy e^V s s ^ t to p a (x) — p% 
is drawn in Fig. |3[ 

The electrostatic potential decays only as exp(— k d x) at large distances 
from the wall, 

$(x) ~ $ as e- KDX (6.1) 

x— >+oo 
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The sign of <£ as depends on the composition {e 7 ,p^} 7= i r .. !ns as well as on 
n D b and A e i. Indeed 

$ as = -— V e*p* {M^Knb, Ad) + [M 7 - M] (x = 0; e 7 ,K D b, A cl )} 

(6.2) 

where, in the limit r\ = K D b = 0, M as is obtained from ( p. 24 ) 



7T 



M„(„ = 0, A ol ) = ^ r ^!_ |__ + (1 + 2Ao]) ln 3 



(1 - A 2 ) ~\ 
+2 { . Clj ln(l-A el ) (6.3) 

Ael J 

while [M 7 — M] (x = 0; e 7 , 77 = 0, A e i) is calculated from ( |3.23 ). 

In an asymmetric plasma or a charge-symmetric plasma with at least 
four components and in an asymmetric state, as far as the density profile of 
one species is concerned, the exp(— 2k d x)/x tail of the screened self-energy 
is always overcome by the effect of — e a $(x), 

P* ym W ~ p^[l-/3e^ ase — ] (6.4) 

x— >+oo L J 

On the contrary, in a symmetric two-component plasma or in a charge- 
symmetric plasma with more than two components and in a symmetric 
state (see definition (|4.4[) ), $(z) = and 



B 



a;— >+oo 



l-/?e 



4.r 



(6.5) 



However, because of the bulk local neutrality, the influence of <J?(x) at 
large distances disappears in the total particle density p a {x) even in the 
case of an asymmetric plasma 



K 2 e -2K D x / e -3hi D x 



z — ' x— >+oo z — ' 1D7T X \ t / 



a 



The total particle density is submitted to an effective repulsion far away 
from the wall. On the contrary, the charge density at large distances is 
determined by $(x), 



H 2 



Te aPa (x) ~ --HQ^e-"** (6.7) 

^— ' x— >+oo 47T 
a 

When approaching the bulk region, the charge density vanishes with a sign 
ruled by the composition of the Coulomb fluid. 
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6.2 Effect of the wall dielectric response 

The three effects which interplay in the density profiles when e w 7^ 1 have 
been discussed in subsections 5.3 and 6.1. Here we consider an asymmet- 
ric two-component plasma with e+ = — 2e_ and we comment briefly the 
corresponding figures for the profiles of particle and charge densities. 

First we turn to the density profiles. When e w = 1 the density at the 
wall (x = b) differs from its bulk value only by a term of order e D = n D f5e 2 /2 
(see Fig. ^:b). When e w < 1 the electrostatic repulsion from the wall makes 
all density profiles at the wall (x = b) vanish exponentially fast when b 
goes to zero (see Fig. f|:a). When e w > 1 the density at x = b increases 
as e w gets larger because of the electrostatic attraction to the wall (see 
Fig. ^}b). Subsequently, the difference p+(x) — p+, which has a double-layer 
structure when e w < 1, exhibits a threefold-layer structure when e w becomes 
sufficiently large. 

The charge density profile C[x) = e + p + {x) + e-p-(x) obeys the same 



evolution when e w varies. If e w = 1, according to (5.6), the condition p+e+ + 
P-e^ > implies that Q HW (x = b) < and enforces the double layer ©ffi 
shown in Fig.|| and discussed after ( |5.7|) . This double layer arises from the 
balance between the electrostatic force associated with and created by 
C (x) itself and the geometric repulsion from the wall due to the deformation 
of screening clouds. When e w < 1, the extra electrostatic repulsion from the 
wall does not destroy the double layer and only enforces the vanishing of 
C(x) at x = when 6 = (see Fig. ^:a). When e w > 1, in the case 
of positive charges, the electrostatic self-attraction to the wall competes 
with the opposite effect of <I>(x) and for high enough values of e w , C(x) 
contains three layers © © (see Fig. |(§b). When e w becomes far larger, 
the electrostatic self-attraction to the wall is so strong that the sign of <£(x) 
at large distances changes and again there appears a double layer © © but 
with signs reversed with respect to the situation when e w = 1 (see Fig. ^:c). 

A Appendix 



In the present appendix we show that the contact theorem ( 4.1 2| ) is satisfied 



by the profile densities found at first order in £ D . First, we consider a plain 
hard wall (A e i = 0) and then we turn to the case of a wall with a dielectric 
response. 

In the case of a plain hard wall, the r.h.s. side of ( [4,12| ) is reduced to 
the kinetic pressure on the wall 

a 

According to the bulk local neutrality ( |2.3j ) the contribution from $(6) to 
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Pa | a 1= o ( x = ^ s zero ; anci , according to ( |5.11| ), 

E ^Ia o1 =o (* = b) = Y<pI - /*E^ e * K£Ia*=o (* = ft ) 

+ 0(pe£ lnE D ) (A.2) 



Since V£ t \ A (x = 6) = —k d /Q (see Section 5), the kinetic pressure 
on the wall does coincide with the value ( }4.13| ) of the bulk pressure at first 
order in e D . 

Now we turn to the case of a wall with a dielectric response. First, the 
last integral in the r.h.s. side of ( 4.12Q can be written as the sum J w — I B 
where 



dxy^p a (x)J Wta (x) 



(A.3) 



and Jw,a{%) is the integral in the r.h.s. of the first BGY equation ( |2.1| ). 
I B and J B , a (x) are defined in a similar way. The value J$ &(x) of J w ,a(x) 
in the mean-field approximation is given in (|2.7|). By use of the definitions 
(1.4) and ( 2.10 ), it can be rewritten as 



.2 d 



-P<te [V£(x) - V scl[ (x)} 



(A.4) 



Therefore 



P 



[ + °° , V- , , d[ejV sel{ ] 



+oo 



dx y^p a (x 



d [ejV; c 
dx 



<x) (A.5) 



The r.h.s. of equality ( |A.5| ) does not involve the bare self-energy V sel f(x) 
and it is calculated by inserting the expression (|l.8| ) of p a (x) in terms of the 
screened self-energy and the electrostatic potential <J>(x). After an integra- 
tion by parts, we notice that, according to ( |3.15| ) and ( p. 16 ), 



dxexp[-0elV^(x)] 



dx 



dx- 



dx 



x [l + 0(e D ]ne D )] 
(A.6) 



Then the local neutrality in the bulk ( |2.3| ) implies that the sum ( |A.5| ) is 
equal to 



E P ° ~ E = & ) + ° (P £ c ln£ i 



(A.7) 
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Second, we calculate I B which is defined in terms of 
JB,a = -(3e a J dx J dy -^{x - x ,y)^2e 1 p 1 {x')h ai {x,x ,y) (Ai 



as in (A.3). In fact, the difference 

h a >y(x,X , y) — ha-y\j\ q (x,x',y) 



(A.9) 



0:7 
is 



gives a vanishing contribution to I B . Indeed, the mean-field value ( |2.5| ) of h, 
is proportional to 4>^ and, according to (|3.2| ) and fl3.4| ), c/>(°) — . 
bounded by a product of functions f(x, x')g(y) where <?(y) is integrable while 
f(x,x') decays exponentially fast in all directions in the plane of variables 
(x,x ! ). As a consequence 



dv 



-^-(x-x',y)^2 e 2 a e 2 p a (x)p^(x') 



a,7 



<f>M(x,x',y)-ct>M ( x ,x',y) 



(A.10) 



is absolutely integrable in the space (x,x',y). Moreover, — (f)^\ A _ n 
is symmetric under exchange of x and x' whereas dv B /dx is antisymmetric 
under the same exchange. Subsequently, the contribution from the difference 
(A.9) to I B is just zero. Moreover, when A e i = the self-energy is bounded 
for all x's ranging from to +00 so that, according to ( 3.17| ), 



Pa{x) - Pa\ Acl=0 (x) = pl {exp [- (3e a V^(x)] - 1 + O exp (e D )} (A.ll) 

and, according to (|3.15| ) and fl3,16|) , p a {x) — p a \^ 1=0 (x) contributes to I B 
by a term of order pe 2 D In well as p-y(x') — p 7 | A =0 (x'). Eventually, 



J b = ^b| Ac1 =o + ° (P £ l ln£ o) 



(A.12) 



Since the Ursell function h ai \ A 1= o ODe y s the BGY hierarchy with v w re- 
placed by v B , a calculation similar to that performed for I w leads to 



BlA rl =0 



b) + 0(pe 2 D lne D ) (A.13) 



According to ( A. 2 ), the r.h.s. of ( A.13| ) coincides with the opposite of the 
term of order e D in the bulk pressure. Therefore ( |A.13| ) leads to 



h = ^ + 0{pe 2 D \ne Dj 



(A.14) 



Eventually, the sum of the terms in the r.h.s. of (4.12) does coincide with 
the value (g]J ) of (3P B . 
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Figure 1: Profile of the electrostatic potential <J?(x) for n D b = 0.1 in the 
limit A e i/3e 2 /6 <C 1. Af~ in ( 3.23 ) does no longer depend on 7 in the limit 
A e i/3e 2 /6 -C 1. $0 = — Ylj P e yPy/ K D- The values of e w are displayed in the 
figure. 




Figure 2: a Figure 2: b 

Figure 2: Density profiles in a two-component plasma near a plain hard wall. 
If the plasma is symmetric (Fig. ^:a) <& HW (x) = and only the geometric 
repulsion from the wall is involved. The curve is the same as in Ref.||]. If 
the plasma is charge asymmetric with e+ = 2|e_| (Fig. ||:b), <fr HW (x = b) < 
= § HW (x = +00) and the competition between the geometric repulsion 
from the wall and the attraction to the wall by <& HW (x) for positive charges 
results into a double layer structure for p+(x) — p\. In Figs ||:a and ^:b 
Ko f3e 2 = 0.01 and n D b = 0.1. 
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Figure 3: a Figure 3: b 

Figure 3: Contribution from the complete screened self-energy e^,V^ elf to 
the profile density. In the case of a symmetric two-component plasma of 
charges e and — e, p+[x) = p-(x) = p(x) and <£(x) = while p(x) = 
p B exp [— Pe 2 V a a ^ { (x)] . In Fig. |3]:a, where e w < 1, both the electrostatic and 
geometric repulsions from the wall make p(x) < p B . Fig. [jjb, where e w > 1, 
displays the competition at short distances between the electrostatic attrac- 
tion to the wall (which gets larger when e w increases) and the geometric 
repulsion from the wall. In Fig. |3|:a, K D (3e 2 = 0.1 and n D b = 10~ 5 whereas 
in Fig. [|:b, n D f3e 2 = 0.01, K D b = 0.1 and the values of e w are written in the 
figure. 
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Figure 4: a Figure 4: b 

Figure 4: Density profiles in an asymmetric two-component plasma (e+ = 
2|e_|) when e w varies. In Fig. §a, K D j3e 2 = 0.1 and k b 6 = 10~ 5 , whereas 
in Fig. ||]:b, K D f5e 2 = 0.01, k d 6 = 0.1 and the values of e w are given in the 
figure. 
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asymmetric 



K D (x-b) 

Figure 5: Double layer of the charge density profile C{x) = e + p + (x) + 
e_p_(j;) for the asymmetric two-component plasma already considered 
in Fig. ||. X} 7 P7 e 7 > in this case and inequality ( p 7 ^ ) can be checked. 
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Figure 6: c 



Figure 6: Charge density profiles in the same two-component plasma as in 
Fig |5[ In Fig. |(|a , K D (3e 2 = 0.1 and K D b = 1CT 5 and in Figs. ^:b and ||:c, 
n D (3e 2 = K D b = 0.1. 
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